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Abstract 

A  computational  procedure  and  some  numerical 
results  of  unsteady  interaction  of  a  helicopter 
rotor  blade  with  a  Lamb-like  vortex  of  finite 
viscous  core  in  subsonic  and  transonic  flows  is 
presented.  The  interaction  considered  here  is  one 
of  the  limiting  cases  of  a  more  complex  interac¬ 
tion  typically  encountered  on  a  helicopter  rotor 
blade.  In  this  limit,  the  interacting  flow  field 
is  considered  to  be  unsteady  but  two-dimensional. 
Accordingly,  unsteady,  two-dimensional,  thin-layer 
Navier-Stokes  equations  are  solved  using  a 
prescribed-vortex  method  (also  called  perturbation 
method)  for  the  cases  of  stationary  and  moving 
rotor  blades  encountering  a  moving  vortex  passing 
the  blades.  The  numerical  results  are  compared 
with  the  recent  experimental  data  of  Caradonna 
et  al.  for  the  latter  case.  The  comparison  shows 
that  for  the  transonic  cases,  the  flow  field  is 
dominated  by  the  presence  of  the  shock  waves,  with 
strong  indications  of  unsteady  time  lags  in  the 
shock-wave  motions  and  shock-wave  strengths,  and 
of  important  three-dimensional  effects.  For 
subcritical-flow  cases,  however,  the  unsteady  lag 
effects  on  the  basic  rotor  blade  are  absent,  and 
three-dimensional  effects  appear  to  be  negligible, 
unlike  the  supercritical  case.  The  subcritical 
calculations  are  in  good  agreement  with  the  exper¬ 
imental  data The  presence  of  three-dimensional 
influence  in  cfte  supercritical-flow  case  is  demon¬ 
strated  by  making  calculations  with  a  small- 
disturbance  code,  After  i •'eluding  in  it  the  neces¬ 
sary  cross-flow  termh.  x'his  modification  to  the 
two-dimensional,  transonic  small-disturbance 
formulation  to  produce  the  three-dimensional 
influence  clearly  demonstrated  the  presence  of 
such  effects,  and  these  numerical  results  are  in 
good  qualitative  agreement  with  the  experimental 
data.  A  similar  calculation  done  for  the 
subcritical-flow  condition  showed  almost  no  three- 
dimensional  Influence,  suggesting  that  the  three- 
dimensional  influence,  if  any,  seems  to  be 
accentuated  in  the  presence  of  shock  waves.  In 
all  these  calculations,  the  details  of  the  vortex 
structure  were  found  to  be  very  important  in 
matching  the  experimental  data. 
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Nomenclature 

A,B,C^ ,C2,D,E  s  coefficients  as  shown  in  Eq.  (10) 
A,B  =  Jacobian  matrices 

an  =  vortex  core  radius 
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=  free-stream  sound  speed 

s  F’(z)/F(x);  see  Appendix 

=  characteristic  length  scale,  chord 
of  the  rotor  blade 

=  lift  coefficient 

:  quarter-chord  pitching-moment 
coefficient 

=  Coefficient  of  pressure 

*  chord  of  the  vortex  generating 
wing 

=  flux  vectors 

=  total  energy  per  unit  volume 

s  flux  vectors 

=  identity  matrix 

=  transformation  Jacobian 

-  Jacobian  matrix  of  viscous  stress 
flux  vector 
=  local  Mach  number 

=  blade  reference  Mach  number 

(Qrg/a^) 

=  blade  tip  Mach  number 
=  free-stream  Mach  number 
=  static  pressure 


=  velocity  Induced  by  the  vortex, 
Iuv  ♦  K 


-  unknown  flow-field  vector 


Euler  solution  of  vortex  In  a 
uniform  free  stream 


angle  between  vortex  axis  and 
rotor-blade  chord  line 


radius  of  rotor  blade 
Reynolds  number 

radial  distance  from  the  vortex 
center 

reference  station  on  the  rotor 
blade,  fraction  of  R 

initial  vortex  position, 

Ixo  + 


viscous  flux  vector 


contravariant  velocity  components 

free-stream  velocity,  forward 
velocity  of  rotor  Tu  ♦  Jv 

velocity  components  in  physical 
plane 

velocity  components  induced  by  the 
vortex 

tangential  velocity  of  vortex 

initial  vortex  location  in  the 
flow  field 

instantaneous  vortex  position 

physical  plane  coordinates 

defines  airfoil  surface 
Yb  1  yb(x-fc) 


angle  of  attack 
vortex  strength 

dimensionless  strength  of  vortex 
(r/U.C) 

ratio  of  specific  heats 

airfoil  thickness  ratio 

inverse  of  blade  aspect  ratio 

explicit  and  implicit  smoothing 
coefficients 

angle  between  the  vortex  velocity 
vector  and  y-axls 
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-  advance  ratio  and  local  advance 
ratio,  respectively 

•  metrics  of  transformation 

=  transformed  plane  coordinates 

=  density 

=  disturbance  potential 
=  azimuth  angle 

=  angular  velocity  of  the  rotor 
blade 

=  reduced  frequency 
=  Indicates  dimensional  value 


Subscripts 


=  free-stream  condition 


1 .  Introduction 

One  of  the  important  problems  of  helicopter 
aerodynamics  and  one  that  was  the  subject  of  many 
recent  experimental 1-11  and  theoretical®-11  studies 
is  the  mechanism  of  rotor  blade  and  vortex  inter¬ 
action.  This  interaction  mechanism  is  a  primary 
source  of  impulsive  noise  generation.  The  tran¬ 
sonic  speeds  of  today's  helicopter  blade  tips  adds 
to  the  complexity  of  the  problem.  The  blade  tips, 
which  trail  the  strong  and  concentrated  tip 
vortices  in  such  a  flow  field,  trace  out  prolate 
cycloidal  paths  in  space,  and  in  the  process 
encounter  a  variety  of  blade-vortex  interac¬ 
tions.  These  interactions  Induce  unsteady  blade 
loading  and  aerodynamic  noise,  with  compressibil¬ 
ity  playing  an  Important  role  in  the  problem. 

Although  the  generic  problem  of  the  blade- 
vortex  interaction  can  be  viewed  in  general  as 
unsteady  and  three-dimensional,  in  one  limit,  when 
the  Intersection  angle  of  the  vortex  with  the 
blade  is  very  small  or  zero,  the  interaction  can 
be  approximated  to  be  two-dimensional  but 
unsteady,  as  shown  schematically  in  Fig.  1.  Most 
of  the  recent  numerical  studies®-11  that  have 
addressed  this  problem  have  been  solved  in  this 
limit.  Transonic  small-disturbance  equations, ®-^ 
full-potential  equations,11  Euler  equations,^-10 
and  thin-layer  Navier-Stokes  equations^  have  all 
been  solved  for  the  problem  of  a  convecting  vortex 
passing  and  interacting  with  the  flow  field  of  a 
stationary  airfoil  in  a  uniform  free  stream. 

Some  of  the  above  methods  place  limitations 
on  the  intensity  of  the  interaction  in  terms  of 


vortex  strength,  of  vortex  location  with  respect 
to  the  airfoil,  and  on  the  free-stream  Mach 
number;  nevertheless,  they  all  seem  to  give  simi¬ 
lar  results.  For  flows  dominated  by  the  strong 
viscous  interaction  and  shock-wave/boundary-layer 
interaction,  the  natural  choice  of  the  equation 
set  that  describes  such  a  flow  completely  would  be 
the  Navier-Stokes  equations.  The  present-day 
numerical  algorithms  for  this  set  of  equations  are 
still  very  expensive  in  terms  of  computing  time, 
although  good  progress  has  been  made  toward 
reducing  the  time  required.  But  the  memory  of  the 
available  computers  is  adequate  to  address  these 
kinds  of  problems  in  the  two-dimensional  limit. 
This  is  particularly  true  with  the  present 
prescribed-vortex,  or  perturbation,  method  which, 
has  been  demonstrated  to  resolve  important  flow 
features  in  a  blade-vortex  Interaction  problem  in 
both  subsonic  and  transonic  flows,  even  in  a  very 
sparse  finite-difference  mesh.^'^ 

A  problem  that  is  more  practical  and  numeri¬ 
cally  more  complicated  than  the  one  mentioned 
above  is  that  of  a  rotating  blade  of  a  helicopter 
rotor  encountering  a  vortex  generated  upstream. 
Such  an  experiment  was  done  recently  in  a  wind 
tunnel  at  the  U.S.  Army  Aeromechanics  Laboratory* 
at  Ames  Research  Center.1  A  schematic  of  the 
experiment  is  shown  in  Fig.  2.  A  vortex  generated 
at  the  tip  of  a  straight  NACA  0015  wing  interacts 
with  the  flow  field  of  a  rotating,  two-bladed 
helicopter  rotor  blade  under  subsonic  and  tran¬ 
sonic  flow  conditions.  For  such  a  rotating  blade, 
unlike  the  stationary  airfoils  considered  in  pre¬ 
vious  studies, 5*11  the  unsteady  time-lag  effects 
that  are  present  even  in  the  absence  of  the  inter¬ 
acting  vortex  are  very  important.12  The  effect  of 
this  unsteady  time  lag  is  to  delay  the  process  of 
adjustment  of  the  flow  corresponding  to  the  appro¬ 
priate  azimuthal  blade  position.  This  has  a  pro¬ 
found  influence  on  the  unsteady  blade  loads.  At 
transonic  conditions  this  will  be  compounded  by 
the  presence  of  strong  shock  waves.  In  the 
results  to  be  presented  later,  it  will  be  demon¬ 
strated  that  these  unsteady  time-lag  effects  are 
Important  only  when  the  flow  is  supercritical; 
they  are  less  Important  if  the  flow  is  subcrlti- 
cal.  In  addition,  the  results  may  also  be  influ¬ 
enced  by  the  three-dimensional  (tip)  effects,  if 
the  chordwise  reference  station  considered  is 
closer  to  the  tip  of  the  blade.  Above  all,  one 
other  very  important  ingredient  of  the  vortex 
Interaction  study  Is  the  detailed  knowledge  of  the 
structure  of  the  interacting  vortex.  Although  an 
analytical  representation  of  the  vortex  is  often 
used,  the  point  to  bear  in  mind  is  that  such  a 
representation  should  have  an  accurate  core 
structure  embedded  in  it. 

In  the  present  investigation,  from  the 
experimental  details  of  Ref.  1,  we  find  that  the 
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rotor-blade  reference  station  (section  AA  in 
Fig.  2)  that  undergoes  parallel  blade-vortex 
interaction  is  inboard  of  the  blade-tip  region. 

If  the  tip  influence  is  neglected  at  this  station, 
the  flow  can  be  approximated  as  two-dimensional, 
but  unsteady.  One  important  difference  of  this 
problem  from  the  earlier  formulation^  is  that 
the  blade  is  rotating  now.  After  implementing 
this  important  feature,  along  with  the  necessary 
changes  in  the  boundary  conditions,  unsteady,  two- 
dimensional,  thin-layer  Navier-Stokes  equations  in 
3trong  conservation-law  form  are  solved  for  the 
interaction  flow  field  of  the  rotating  blade, 
using  an  approximately  factored,  implicit  finite- 
difference  numerical  algorithm  written  in  delta- 
form.  An  analytical  representation  of  the 

measured  vortex  structure  was  used  in  the 
computations. 

Although  the  initial  finite-difference  grids 
were  generated  by  an  algebraic  method,11*  an 
adaptive-grid  procedure'®’1®  was  used  throughout 
so  as  to  resolve  the  important  flow  features, 
including  shock  waves,  vortex/shock- induced  sepa¬ 
ration  of  the  boundary  layers,  if  any,  and  the 
vortex  structure  itself. 

In  this  paper  the  governing  equations  and 
numerical  formulations  are  discussed,  and  the 
numerical  results  and  comparisons  with 
experimental  data  are  presented. 


2.  Governing  Equations  and  Solution  Procedure 

The  interaction  flow  field  is  solved  bv  the 
prescribed-vortex,  or  perturbation,  method.  '  The 
essence  of  the  method  is  that  each  of  the  depen¬ 
dent  flow  variables  is  split  into  a  prescribed 
part,  which  is  simply  the  vortical  disturbance, 
and  a  remaining  part,  which  is  obtained  from  the 
solution  of  the  governing  equation  set.  Even 
though  the  governing  equations  are  nonlinear  and 
independent  solutions  are  not  superposable,  we  may 
still  decompose  the  dependent  variables  as 


q  =  %  *  (q  -  V 


ou  and  qo  =  ou 
pv  pv 


Here  q  is  the  unknown  flow-field  vector,  and  the 
vector  q  represents  the  solution  of  Euler  equa¬ 
tions  for  the  vortical  disturbance  convecting  in  a 
uniform  free  stream.  Previous  studies^*^  have 
demonstrated  that  the  perturbation  method  can  in 
fact  resolve  the  flows  with  concentrated  vortices 
better,  even  with  a  coarser  grid,  than  can  a  non¬ 
perturbation  method.1®  For  example,  a  representa¬ 
tive  plot  of  the  variation  of  lift  coefficient  as 
a  function  of  the  vortex  position  is  reproduced 
from  Refs.  8  In  Fig.  3.  This  calculation  was  done 


using  Guler  equations  for  the  case  of  a  oonvecting 
vortex  encountering  a  nonlifting,  stationary  NACA 
64A006  airfoil  in  transonic  flow,  and  using  the 
sane  grid  topology  with  both  perturbation  and 
nonperturbation  (or  conventional)  methods.  This 
figure  clearly  demonstrates  that  with  the  con¬ 
ventional  nethod  the  accuracy  of  the  solution 
increases  with  the  fineness  of  the  finite- 
difference  mesh,  whereas  the  perturbation  method 
produces  a  much  better  solution  even  in  a  coarser 
mesh. 


dependent  variables  p,  p,  e,  u,  and  v  by  p^, 
tpa,  a„,  respectively. 

The  pressure,  density,  and  velocity  compo¬ 
nents  are  related  to  the  energy  per  unit  volume  by 
the  equation  of  state,  which  is  written  for  a 
perfect  gas  as 


-b*-(2LF!) 


The  governing  partial  differential  equations 
are  the  unsteady,  thin-layer  Navier-Stokes  equa¬ 
tions.1^  The  equations  are  written  in  nondimen- 
sional,  strong  conservation-law  form  for  a  perfect 
gas  using  the  generalized  Independent  coordinate 
system  of 

5  s  5(x,y,t) 

n  =  n(x,y,t)  (2) 

t  =  i(t) 

and  in  the  perturbation  form^*^  as 

3t(q  -  qo>  *  3e(E  -  Eq)  ♦  -  Fq)  =  Re  3^ 


where  '  denotes  quantities  scaled  by  the 
Jacobian,  for  example 


This  equation  of  state,  along  with  the  mass  and 
momentum  equations  given  by  Eq.  (3),  complete  the 
equation  set  to  be  solved. 


*  -1 
q  =  J  q  ; 


J  =  Vy  *  V*  1  1/(Vn  *  V?5 

is  the  transformation  Jacobian.  The  flux 
vectors  E,  E  ,  F,  F  ,  and  S  are  described  in 

O  O  A 

detail  In  Refs.  7-9.  The  viscous  flux  vector  S 
is  written  in  the  context  of  a  thin-layer  model ^ 
and  hence  is  valid  for  high-Reynolds-number  turbu¬ 
lent  flows.  The  turbulent  eddy  viscosity  is  com¬ 
puted  U3lng  a  two-layer,  algebraic  eddy-viscosity 
model.20 

The  generalized  coordinate  system  of  t,n,x 
allows  the  boundary  surfaces  in  the  physical  plane 
to  be  mapped  onto  rectangular  surfaces  in  the 
computational  plane.  This  feature  simplifies  the 
procedure  of  grid-point  clustering  in  the  flow 
regions  that  experience  rapid  change  in  the  flow- 
field  gradients. 

The  primitive  variables  that  make  up  the 
governing  equations,  Eq.  (3),  are  density  p,  the 
two  mass  fluxes  pu  and  pv  in  the  two  coordinate 
directions  x  and  y  (where  x  is  the  streamwise 
direction  and  y  is  normal  to  It),  and  the  total 
energy  per  unit  volume  e.  All  length  scales  are 
normalized  by  the  chord  (assumed  to  be  unity)  of 
the  rotor  blade  at  the  reference  station,  and  the 


In  the  present  formulation,  the  difference 
between  the  problems  of  the  stationary  blade  and 
moving  blade  show  up  in  the  metric  terms  involving 
time  derivatives,  namely,  and  n..  Following 
the  formulation  given  by  Isom21  and  Caradonna  and 
Isom22  for  unsteady  flow  over  helicopter  rotor 
blades,  the  effective  local  Mach  number  at  the 
reference  station  can  be  written  as 

M.  =  M_(  1  ♦  ii’  sin  wt)  (6) 

l  P 

where  Mp  is  the  rotational  Mach  number,  w'  is 
the  local  advance  ratio,  and  u  is  the  effective 
reduced  frequency.  This  formulation  enables  the 
flow-field  solution  of  a  rotating  blade  with  an 
oncoming  free  stream  to  be  solved  in  the  two- 
dimensional  limit  as  a  blade  moving  with  an 
x-velocity  of  Mra_p'  sin  ut  in  a  flow  of  Mach 
number  Mr.  Note  that  both  and  nfc  are  non¬ 
zero  for  this  case,  whereas,  for  a  stationary 
blade  with  a  fixed  grid,  both  and  are 
zero. 

The  boundary  conditions  are  applied  explic¬ 
itly.  Since  the  grid  extends  20  chord-lengths  in 
all  directions  from  the  surface  of  the  blade, 
free-stream  conditions  are  specified  at  the  outer 
boundary,  and  simple  extrapolation  is  used  for  p, 
pu,  and  pv  at  the  outflow  boundary.  For 
supersonic  flow,  the  total  energy  e  is  also 
extrapolated,  but  for  subsonic  flow,  the  pressure 
is  held  constant  at  the  free-stream  value,  and 
e  is  obtained  from  Eq.  (5).  To  ensure  continuity 
across  the  wake  cut,  the  flow  variables  are 
linearly  extrapolated  to  obtain  the  values  along 
the  cut. 

Along  the  body  surface  n(x,y,t)  =  0,  the 
no-slip  condition  for  viscous  flow  without  suction 
or  injection  is  given  by  setting  U  =  0  and 
V  =  0.  The  pressure  along  the  body  surface  is 
obtained  by  solving  the  normal  momentum  equation, 
and  the  density  at  the  surface  is  obtained  by 
extrapolation  from  the  grid  interior.  The  total 
energy  e  is  calculated  from  the  known  pressure 
and  density  at  the  surface.  The  boundary  condi¬ 
tions  are  of  low  order,  and  hence  require  that  the 
grid  be  clustered  and  normal  at  the  body  sur¬ 
face.  The  interacting  vortex  is  initialized,  as 
in  Ref.  9. 
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Surface-conforming  grids  are  needed  to 
simplify  the  application  of  the  body-boundary 
condition  procedure.  In  this  study,  an  adaptive 
gridding  procedure  of  Nakahashi  and  Deiwert  J’ 
is  used  to  resolve  flow  features  and  to  Improve 
the  accuracy  of  the  numerical  method.  Briefly, 
the  method  uses  tension  and  torsion  spring  analo¬ 
gies.  The  tension  spring,  which  connects  the 
adjacent  grid  points  to  each  other,  controls  grid 
spacings  so  that  clustering  is  obtained  in  regions 
containing  shock  waves  and  shear  layers.  On  the 
other  hand,  the  torsion  spring,  which  is  attached 
to  each  grid  node,  controls  inclinations  (angles) 
of  coordinate  lines  and  prevents  excessive  grid 
skewness.  The  mesh  can  be  made  nearly  orthogonal 
at  the  surface.  A  marching  procedure  is  used  that 
results  in  a  simple  tridiagonal  system  of  equa¬ 
tions  at  each  coordinate  line  to  determine  the 
grid-point  distribution.  Multidirectional  grid 
adaptation  is  achieved  by  successive  application 
in  each  direction.  For  the  compressible  flow 
fields  considered  in  this  study,  the  density  gra¬ 
dient  was  found  to  be  the  best  choice  to  drive  the 
adaption  in  the  x-direction;  the  Mach  number 
gradient  was  the  best  choice  of  driver  in  the 
normal  direction.  In  actual  practice,  given  a 
baseline  grid,  the  above  procedure  will  modify  the 
grid  at  specified  Intervals  to  resolve  the  flow 
satisfactorily.  In  the  present  study,  the  grid 
was  adapted  at  every  two-marching-step  inter¬ 
vals.  This  increased  the  computational  time  by 
about  50%.  The  baseline  grid  used  was  a  surface- 
conforming  C-grid  generated  by  an  algebraic 


method  of  Pulliam  et  al.11*  and  had  221  points 


around  the  airfoil  and  67  points  in  the  normal 
direction.  The  grid  boundary  was  chosen  to  be  at 
20  chords  in  all  directions. 


An  implicit,  spatially  factored  numerical 
algorithm  with  Euler-implicit  time-differencing 13 
is  used  to  solve  the  perturbation-form  of 
Eq.  (3).  This  algorithm  is  written  in  delta  form 
as 
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where  A,  B,  and  M  are  the  Jacobian  matrices 
detailed  in  Ref.  19;  I  is  the  identity  matrix; 
4.,  fi_  are  the  spatial  central  difference  opera- 
4  and  V  are  the  forward  and  backward 


dr 


n  . 
rs;  and 


difference  operators.  For  convenience 
45  =  1  =  fin  is  assumed.  The  time  index  is 

denoted  by  h,  and  qn  =  qn(nfit),  fiqn  =  qn+1  -  qn 
and  cT  and  cp  are  the  implicit  and  explicit 
ng  coefficients,  respectively.  Second- 
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order  implicit  and  fourth-order  explicit  numerical 
dissipation  terms  are  added  to  the  numerical 
scheme  to  improve  the  nonlinear  stability  limits 
posed  by  thefine  mesh.2^  Even  so,  the  nondimen- 
sional  time-steps  generally  were  restricted  by  the 
stability  constraints  to  the  order  of  0.05°  of  the 
azimuthal  travel  of  the  blade  motion  or  0.005455 
of  the  chord  travel . 


The  numerical  scheme  is  first-order  accurate 
in  time  and  second-order  accurate  in  space.  Fur¬ 
ther,  in  writing  Eq.  (7),  it  is  assumed  that 


A  a  A  and  that  B  a  B,  where 


o 

and 
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Central  differencing  is  used  throughout  the 
solution  domain,  except  in  regions  of  supersonic 
flow  before  shock  waves,  where  upwind  differencing 
is  used.  The  transformation  metrics  are  not  known 
analytically  and  are  computed  numerically  by  cen¬ 
tral  differencing  (second-order)  at  the  interior 
points  and  by  three-point,  one-sided  differencing 
at  the  boundaries. 


3.  Results  for  a  Stationary  Rotor  Airfoil 


In  this  section,  numerical  results  are  pre¬ 
sented  for  a  moving  vortex  encountering  a  fixed 
rotor  blade  under  transonic  conditions.  All  cal¬ 
culations  are  done  for  NACA  0012  airfoil  in  a 
uniform  free  stream  of  M^  =  0.8  at  a  =  0°,  and 
assuming  a  turbulent  boundary  layer.  As  mentioned 
earlier,  an  adaptive  gridding  procedure  is  used  in 
all  results  presented  here. 


A  Lamb-like  vortex  with  a  finite  viscous  core 
(aQ  =  0.05)  and  a  cylindrical  velocity  distribu¬ 
tion  given  by^-^ 


vjr) 
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(8) 


was  chosen  to  interact  with  the  flow  field  of  a 
stationary,  nonlifting  rotor  airfoil  in  a  uniform 
free  stream  of  Mach  number  M,,.  The  strength  of 
the  vortex  and  its  location  with  respect  to  the 

airfoil  were  chosen  to  be  r  =  0.2  and 
yv  =  -0.26,  respectively.  The  interaction  flow 
field  was  computed  in  the  same  manner  as  outlined 
in  earlier  studies. ^  Figure  4  shows  plots  of 
Instantaneous  surface-pressure  distributions,  the 
local  grid  arrangement,  the  Mach  number  contours, 
and  the  pressure  contours  for  different 
x-locations  of  vortex  positions  as  the  vortex 
passes  by  the  airfoil.  The  passing  vortex  induces 
on  the  airfoil  a  continually  changing  effective 
angle  of  attack.  Because  of  the  sense  of  rota¬ 
tion,  it  induces  a  downwash  initially  when  it  is 
upstream  of  the  leading  edge,  changing  to  upwash 
as  it  passes  behind  the  airfoil.  This  Induces  a 
continuous  change  in  the  blade-loading  pattern. 

It  should  be  noted  here  that  the  initial  lift  on 
the  blade  is  zero  and  that  any  lift  generated 
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during  the  interaction  is  induced  by  the  vortex. 

As  observed  before, the  maximum  Influence  of 
the  vortex  on  the  airfoil  flow  field  seems  to 
occur  when  the  vortex  is  within  one  chord  of  the 
airfoil. 

Previous  calculations  done  for  the  same  air¬ 
foil  under  identical  conditions2^  used  a  fixed- 
grid  topology,  with  221  «  67  grid  points;  in  that 
study,  shocks  were  not  well  resolved.  In  the 
present  case,  the  adaptive-grid  topology,  shown  in 
Fig.  4,  also  uses  221  «  67  grid  points,  but  is 
clearly  able  to  resolve  all  aspects  of  the  flow. 
Since  the  grid  is  adapted  in  both  the  x-  and 
y-directions,  it  clearly  resolves  shock  waves,  the 
interacting  vortex,  and  the  viscous  layer  at  the 
body  surface.  Large-scale  vortex  and  shock- 
induced  boundary- layer  separation  was  not 
observed,  although  the  formation  of  a  small  sepa¬ 
ration  bubble  and  sufficient  thickening  of  the 
boundary  layer  were  seen  behind  shock  wave  for 
conditions  when  the  vortex  influence  was  maximum. 

As  the  vortex  passes  the  airfoil,  it  encoun¬ 
ters  the  shock  wave  sitting  on  the  surface.  The 
vortex  actually  splits  the  shock  wave  into  a 
triple  shock  wave,  as  is  clearly  seen  from  the 
local  grid  structure,  Mach  contours,  and  pressure 
contours  of  Fig.  4  for  a  vortex  position  past  the 
midchord  of  the  airfoil.  This  feature  was  con¬ 
firmed  by  making  a  Euler  calculation  of  the  same 
flow,  by  using  a  fixed-grid  topology,  and  by 
heavily  clustering  the  grid  in  the  region  of 
Interest.  This  apparent  unsteady  effect  seems  to 
fade  away  as  the  vortex  passes  several  chords 
downstream  of  the  trailing  edge  of  the  airfoil. 
Figure  5  shows  a  plot  of  Instantaneous  lift  and 
pitching-moment  coefficients  as  a  function  of 
the  x-vortex  position  during  the  interaction 
process.  The  lift  coefficient  is  initially  nega¬ 
tive  (because  of  the  sense  of  rotation  of  vortex), 
reaches  a  negative  maximum  for  a  vortex  position 
slightly  upstream  of  the  leading  edge  of  the  air¬ 
foil,  increases  to  near  zero  value  at  xy  *  0.7  C, 
and  stays  near  that  value  for  the  rest  of  the 
interaction  period.  Pitching-moment  changes  are 
maximum  when  the  vortex  is  within  1  C  of  the 
airfoil. 


4.  Results  for  a  Rotating  Blade 

In  this  section,  computational  results  are 
presented  corresponding  to  two  experimental  condi¬ 
tions  of  Caradonna  et  al.1  One  condition  each  of 
subcritlcal  and  supercritical  flows  with  and  with¬ 
out  vortex  encounter  will  be  discussed.  The  cor¬ 
responding  tip  Mach  number  of  the  two  cases  is  0.6 
and  0.8,  respectively,  with  an  advance  ratio 
u  -  0.2  for  both. 

The  schematic  of  the  experimental  arrangement 
is  shown  in  Fig.  2.  The  experiments  were  per¬ 
formed  in  a  wind  tunnel  where  a  NACA  0015  wing 
generated  a  tip  vortex  upstream  of  a  two-bladed, 
model  helicopter  rotor.  The  interacting  vortex, 


of  core  radius  25  mm,  passed  the  reference  rotor 
spanwlse  station  at  r0  =  0.893  and  at  a  distance 
of  yv  =  -0.4  C.  The  chord  of  the  rotor  b'.  ide  was 
152.4  mm,  and  the  diameter  of  the  rotor  was 
2.134  m. 

Structure  of  the  Experimental  Vortex 

Comparisons  of  numerical  results  and  the 
above  experimental  data  will  be  meaningful  only 
if  we  consider  the  correct  structure  of  the 
experimental  vortex.  The  information  about  the 
vortex  core  size  and  the  strength  alone  is  not 
sufficient  to  describe  the  correct  vortex  struc¬ 
ture.  Since  the  vortex  peak  velocity  determines 
to  a  large  extent  the  induced  effects  caused  by 
the  vortex  interaction,  an  analytical  representa¬ 
tion  of  the  cylindrical  velocity  field,  such  as 
that  of  Eq.  (8),  based  Just  on  the  vortex 

strength  r  and  the  core  radius  aQ  may,  in 
fact,  give  a  totally  different  vortex  structure. 
Comparison  of  the  cylindrical  velocity  distribu¬ 
tion  of  the  measured  data  of  Takahashi  and 
McAlister2^  and  of  Orloff  and  Grant2**  with  that 
calculated  from  Eq.  (8)  using  aQ  =  25  mm  and 

r  =  0.31  reported  by  Caradonna  et  al.1  showed  a 
40%  greater  peak  velocity  and  a  different  distri¬ 
bution,  although  it  had  the  inviscid  vortex,  1/r, 
behavior  well  outside  the  viscous  core. 

To  circumvent  this  problem,  we  have  put 
together  all  the  available  NACA  0015  data25'26 
that  were  obtained  under  similar  conditions  into 
one  plot  (Fig.  6)  of  cylindrical  velocity  distri¬ 
bution  as  a  function  of  the  distance  from  the 
vortex  center.  Based  on  some  of  the  arguments 
presented  by  Iversen2^  in  correlating  a  large  mass 
of  data  (particularly  vg  on  turbulent  trail¬ 
ing  vortex  decay,  the  expS-Mmental  data  of  Fig.  6 
could  be  represented  by  an  analytical  expression, 
similar  to  Scully's  viscous  vortex,2”  by 


The  core  radius  aQ  was  taken  to  be  25  mm,  as 
reported  by  Caradonna  et  al.1  and  in  agreement 
with  Refs.  25-26.  With  this  aD,  the  value  of 

r  was  chosen  to  match  the  data  of  Orloff  and 
Grant2”  for  vg  .  This  analytical  fit  to  the 

experimental  da¥txis  also  shown  in  Fig.  6  along 
with  1/r  behavior  of  an  inviscid  vortex.  Based 
on  this  analytical  description,  the  strength  of 
the  vortex  was  found  to  be  only  80%  of  the  value 
quoted  by  Caradonna  et  al . 1 

Subcritlcal  Case 

This  condition  corresponds  to  Hj  :  0.6 
and  u  s  0.2.  For  the  reference  station  at  89.3% 
of  the  rotor  blade  radius,  these  conditions 
translate  to  Mr  =  0.536  and  u'  =  0.223. 
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Reaults  for  Rotor  Alone.  First  we  consider  a 
nonlifting  rotor  in  forward  flight  in  the  absence 
of  vortex  Interaction.  The  objective  of  such  an 
exercise  is  to  determine  the  importance  of  three- 
dimensional,  unsteady  time-lag  effects  as  the 
rotor  sweeps  in  azimuth,  say  from  0°-l80°.  As  the 
blade  rotates  from  0°-l80°,  its  effective  Mach 
number  increases  in  the  first  quadrant,  reaching  a 
maximum  at  the  i|>  =  90  +  A*  position;  it  again 
decreases  to  the  average  vaiue  at  *  =  180°.  In 
going  through  this  change  in  local  Mach  number,  if 
the  flow  on  the  rotor  at  *  =  180°  is  nearly  the 
same  as  that  at  *  *  90  -  &i>,  the  the  rotating 
blade  behaves  as  if  it  were  quasi-steady  and 
quasi-two-dimensional.  If,  on  the  other  hand,  the 
two  flows  are  different,  then  the  associated 
three-dimensional  time-lag  effects  will  have  an 
influence  on  the  vortex-blade  interaction  flow 
field  that  develops  in  the  vicinity  of  the 
s  180°  position. 

Figure  7  shows  a  sequence  of  instantaneous 
surface-pressure  results  at  several  azimuthal 
positions  of  the  blade.  For  this  flow,  the  ini¬ 
tial  local  Mach  number  for  4  =  0°  is  0.536  and 
increases  to  a  maximum  (based  on  the  local  advance 
ratio)  at  *  =  90°  and  again  decreases  to  0.536 
at  i|i  =  180°.  Examination  of  these  pressure  dis¬ 
tributions  indicates  that  for  this  subcritical 
flow  condition,  the  unsteady  time-lag  effects  are 
negligible. 


Supercritical  Case 

This  condition  corresponds  to  a  tip  Mach 
number  of  0.8,  with  an  advance  ratio  of  0.2,  and 
the  blade  reference  station  rB  =  0.893,  as 
before.  This  translates  to  a  reference  Mach 
number  of  0.714  and  a  local  advance  ratio  u'  of 
0.223. 

Results  for  Rotor  Alone.  We  first  present 
results  for  the  rotor-alone  case  (i.e.,  in  the 
absence  of  vortex  interaction).  Figure  9  shows 
plots  of  instantaneous  surface-pressure  distribu¬ 
tions  at  four  rotor-blade  azimuthal  positions. 

For  this  case,  the  initial  lift  on  the  rotor  blade 
(at  *  =  0°)  is  zero,  and  the  flow  is  subcritical; 
but  as  the  blade  rotates,  its  local  effective  Mach 
number  increases,  reaching  a  maximum  at  4>  =  90° 
and  decreasing  in  the  second  quadrant.  As  seen  in 
these  pressure-distribution  plots,  the  shock  wave 
that  develops  as  the  blade  rotates  gets  stronger 
and  moves  toward  the  trailing  edge  of  the  air¬ 
foil.  Even  though  the  effective  local  Mach  number 
reaches  a  maximum  for  $  -  90° ,  the  shock  wave 
continues  to  get  stronger  and  move  toward  the 
trailing  edge  in  the  second  quadrant  before  it 
starts  moving  upstream  toward  the  leading  edge  of 
the  blade.  This  unsteady  time  lag  in  the  growth 
and  decay  of  the  shock  wave  persists,  even  when 
the  rotor  blade  has  passed  the  180°  azimuthal 
position  into  the  third  quadrant. 


Vortex  Encounter  with  a  Rotating  Blade.  To 
compute  the  vortex  interaction  flow  field,  as 
before,  the  vortex  was  initialized  at  the 
41  -  0°  position  of  the  rotor  blade  (the  corre¬ 
sponding  xv  s  -19.638)  as  in  Refs.  7-9.  It  then 
convects  with  the  flow  at  the  effective  "free- 
stream  velocity,"  which  i3  Mpa_( 1  ♦  y'  sin  *)  in 
this  case.  Typical  instantaneous  surface-pressure 
distributions  are  shown  in  Fig.  8  corresponding  to 
several  vortex  positions  as  the  vortex  convects 
past  the  rotor  blade.  Since  the  rotor  is  a  non¬ 
lifting  one,  the  lift  is  initially  zero,  and,  as 
the  vortex  approaches  the  blade,  it  Induces  a 
downwash  and  hence  a  negative  lift.  This  continu¬ 
ously  increases  and  reverses  in  sign,  becoming 
positive  lift  as  the  vortex  passes  the  airfoil. 
Shown  in  this  figure  are  data  from  the  experiments 
of  Caradonna  et  al . 1  The  comparison  of  numerical 
results  and  experimental  data  show  good  agreement, 
both  qualitatively  and  quantitatively;  the  peak 
pressures  on  the  nonvortex  side  of  the  blade  seem 
to  be  underpredicted.  It  should  be  pointed  out 
here  that  in  the  experiments,  the  rotor  model  had 
pressure  taps  on  only  one  surface.  To  get  the 
pressure  distribution  on  both  the  surfaces,  the 
model  was  simply  Inverted  and  a  second,  almost 
identical,  experiment  was  run.  Nevertheless,  it 
is  gratifying  to  see  such  a  good  agreement  with 
experiments.  So  far  we  have  neglected  the  three- 
dimensional  effects,  if  any  exist.  Even  if  three- 
dimensional  effects  are  present,  because  of  the 
proximity  of  the  reference  station  to  the  blade 
tip  region,  their  influence  at  this  flow  condition 
appears  minimal. 


Comparison  of  experimental  data2^  with 
numerical  results  in  Fig.  9  shows  relatively  good 
agreement  as  far  as  the  pressure  levels  are  con¬ 
cerned  for  all  azimuthal  positions  of  the  blade  up 
to  *  *  150°,  but  the  location  of  the  shock  wave 
is  too  far  aft  in  the  second  quadrant.  In  the 
experiments,  the  shock  wave  seems  to  have  col¬ 
lapsed  between  ii  -  150°  and  <i  s  180°.  But  the 
numerical  results  show  the  survival  of  a  strong 
shock  wave  even  at  (i  =  180°,  as  discussed 
below.  This  suggests  the  presence  of  possible 
three-dimensional  effects  besides  the  unsteady 
shock-wave  lag  effects. 

Some  insights  into  possible  three-dimensional 
influences  can  be  obtained  from  the  transonic 
small-disturbance  equation  for  a  thin,  high- 
aspect-ratio  rotating  blade21: 

fl*tt  ♦  B*xt  1  h  [C1*x  ♦  C2*xJ 

*  *yy  +  D*xz  ♦  E*zz  (10) 

As  discussed  in  the  Appendix,  this  equation 
reveals  by  inspection  that  three-dimensional 
effects  influence  the  solution  through  the  two 
terms  D«xz  and  E*zz.  where 

,2/3 


and 


D  =  -2Mpey'  cos  t(z  ♦  u '  sin  t)/« 


E  .  e2/«2/3 


Of  these  two  terms,  the  coefficient  G  in  the 
latter  is  independent  of  Mach  number,  and  E*zz 
is  undoubtedly  important  in  the  immediate  vicinity 
of  the  tip  of  the  rotor  blade.  On  the  other  hand, 
the  former  term  0*xz,  is  clearly  dependent  on 
Mach  number,  and  the  coefficient  D  is  also  a 
strong  function  of  the  blade  aspect  ratio,  the 
azimuth  of  the  rotor  blade,  the  spanwise  station 
along  the  blade,  and  the  local  advance  ratio, 
y'  =  U  /ar„,  of  the  rotor.  Therefore,  it  seems 
likely  "’that  this  term  might  play  the  dominant 
three-dimensional  role  inboard  of  the  tip,  espe¬ 
cially  under  transonic  conditions. 

A  qualitative  estimate  of  the  influence  of 
this  cross-flow  term  has  been  deduced  from  two- 
dimensional  calculations  in  which  D*xz  was  added 
ad  hoc  a3  an  inhomogeneous  "source"  term  to  the 
basic  small-disturbance  code  ATRAN2  (Ref.  30), 
as  discussed  In  the  Appendix.  Namely, 

•„<*,y,*,t>  -  F'(z)«x(x,y,t)  was  evaluated  by 
assigning  an  arbitrary  constant  to  F'/F  at  a 
fixed  z  :  rB,  and  by  evaluating  «x  at  the  pre¬ 
vious  time-step.  It  was  found  that  negative 
values  of  F'/F  tended  to  strengthen  the  shock 
wave  for  0  <  *  <  90  and  to  weaken  it  for 
90  <  *  <  180;  positive  values  of  F'/F  had  the 
opposite  effect.  The  sensitivity  of  the  calcu¬ 
lated  results  to  the  value  of  F'/F  is  illus¬ 
trated  in  Fig.  10,  at  *  *  150"  and  180’  for  the 
example  considered  in  Fig.  9-  Figure  10  shows, 
especially  at  the  180°  azimuthal  position,  a 
drastic  alteration  of  the  shock-wave  strength  and 
position  with  increases  in  the  magnitude  of  the 
cross-r’ow  term. 

Figure  11  shows  a  comparison  of  the  experi¬ 
mental  pressure  distributions,  the  Navler-Stokes 
calculations,  and  the  transonic  small-disturbance 
calculations  with  and  without  the  above  three- 
dimensional  correction.  The  two-dimensional 
ATRAN2  and  thin-layer  Navler-Stokes  results  are  in 
good  agreement,  and  both  calculations  show  the 
persistence  of  the  strong  shock  wave  for 
i>  =  180°.  On  the  other  hand,  the  three- 
dimensional  correction  modifies  the  ATRAN2  results 
into  closer  agreement  with  the  experimental  data, 
and  the  calculated  shock  wave  at  *  =  180°  has 
almost  vanished. 

It  must  be  emphasized  that  the  above  analysis 
is  only  intended  to  demonstrate  the  potential 
importance  of  cross-flow  effects  in  supercritical 
flow,  and  not  to  predict  accurately  such  an 
effect.  Although  it  seems  logical  that  three- 
dimensional  effects  should  depend  on  Mp,  u ' ,  V, 
and  e  in  the  manner  suggested  by  the  coeffi¬ 
cient  E  in  Eq.  (10),  the  separation  of  vari¬ 
ables  *(x,y,z,t)  =  F(z)*(x,y,t)  is  only  an 
assumption,  and  the  choice  of  F’/F  is  purely 
arbitrary.  Nevertheless,  the  results  in 
Figs.  9-11  clearly  demonstrate  that  the  rotor- 
alone  supercritical  case  under  consideration  is 
neither  quasl-two-dimensional  nor  quasi-steady. 
This  stands  in  sharp  contrast  with  the  previous 
subcritlcal  case,  Figs.  7  and  8,  for  which  the 


two-dimensional  numerical  calculations  are  in 
excellent  agreement  with  the  model  rotor  data,  and 
for  which  the  significant  unsteady  effects  are 
only  a  result  of  the  airfoil-vortex  interaction. 

Vortex  Interaction  with  a  Rotor  Blade 
Section. A  previous  attempt’  to  compute  the 
blade-vortex  interaction  flow  field,  for  the 
supercritical  flow  case  without  considering  the 
unsteady  time-lag  effects  produced  unsatisfactory 
results  in  terms  of  agreement  with  experiments. 
Typical  results,  reproduced  from  Ref.  9,  are  shown 
in  Fig.  12.  First  indications  of  poor  agreement 
with  experiments  are  the  absences  of  strong  shock 
waves  which  is  a  consequence  of  neglecting 
unsteady  lag  effects  of  the  rotating  blade.  The 
results  presented  in  Figs.  9-11  show  not  only  the 
Importance  of  unsteady  time-lag  effects,  but  also 
the  three-dimensional  influence  for  this  flow. 

Based  on  the  arguments  presented  above,  the 
blade-vortex  interaction  flow  field  at  this  super¬ 
critical  condition  is  calculated,  using  the  two- 
dimensional,  transonic  small-disturbance  equations 
with  the  three-dimensional  correction.  As  shown 
in  Fig.  10,  the  calculations  are  sensitive  to  the 
value  of  the  three-dimensional  correction  factor 
B1.  From  the  results  presented  in  Figs.  9-11,  a 
value  of  B1  s  0.2  gave  the  best  match  with 
experimental  data2®  for  the  rotor  azimuthal  posi¬ 
tion  of  180°.  This  value  for  B1  was  used  in  the 
blade-vortex  interaction  calculations. 

Figure  13  shows  the  results  of  these  calcula¬ 
tions  compared  with  the  experimental  data.1  The 
results  are  in  good  qualitative  agreement. 

Although  the  calculations  show  good  quantitative 
agreement  with  the  experiments  at  the 
xv  =  -0.35  station,  the  agreement  seems  to 
degrade  progressively  as  the  vortex  passes  down¬ 
stream  of  the  blade.  It  should  be  noted  that  the 
flow  on  the  blade  i3  3uberitical  (with  zero  lift) 
in  the  absence  of  the  vortex  interaction  for 
*  s  180"  (see  Fig.  11).  With  the  vortex  interac¬ 
tion,  the  flow  on  the  blade  becomes  supercritical, 
with  strong  shock  waves  as  the  interaction  effects 
peak.  As  the  vortex  passes  downstream  of  the 
blade,  these  supercritical  conditions  progres¬ 
sively  change  to  subcritlcal  conditions.  Also, 
because  of  the  vortex  interaction,  the  blade 
develops  lift  which  is  initially  negative  (because 
of  the  sense  of  rotation  of  the  vortex),  becoming 
positive  as  the  vortex  passes  the  blade. 

In  view  of  the  strong  viscous  interaction 
Inherent  in  this  present  problem,  it  is  not  sur¬ 
prising  to  see  less-than-perfect  agreement  between 
calculations  and  experimental  data.  Probably  this 
is  a  limitation  of  this  present  method.  Neverthe¬ 
less,  it  has  demonstrated  the  importance  of  the 
three-dimensional  effects  in  the  presence  of  shock 
waves.  An  additional  important  observation  indi¬ 
cates  that  if  the  calculated  results  are  lagged  by 
about  2.5°  of  azimuthal  position  (or  about  a 
quarter  chord  of  travel),  the  calculations  seem 
to  show  better  quantitative  agreement  with 


experiments  for  all  x-positions  of  the  vortex. 

--Typical  pressure  distributions,  replotted  with 
this  phase  shift  of  a  quarter  chord,  are  shown  in 
Fig.  14. 

5.  Concluding  Remarks 

This  paper  presents  a  computational  procedure 
for  calculating  the  interaction  flow  field  of  a 
passing  vortex  with  a  helicopter  rotor  blade  in 
forward  flight,  as  well  as  numerical  results  for 
subsonic  and  transonic  flow  conditions.  The 
interaction  of  the  vortex  considered  here  is  one 
of  the  limiting  cases  of  a  more  complex  interac¬ 
tion  typically  encountered  on  a  helicopter  rotor 
blade,  and  it  corresponds  to  the  parallel  blade- 
vortex  interaction  experimental  conditions  of 
Caradonna  et  al.'  In  this  limit,  the  interaction 
flow  field  can  be  approximated  to  be  two- 
dimensional  and  unsteady.  The  numerical  scheme 
involves  the  solution  of  unsteady,  two- 
dimensional,  thin-layer  Navier-Stokes  equations 
implicitly  using  the  prescribed- vortex,  or  pertur¬ 
bation,  approach  presented  in  Refs.  7-9. 


The  influence  of  three-dimensional  effect  was 
further  examined  by  using  a  modification  of  tb= 
ATRAN2  code  (two-dimensional,  transonic,  small- 
disturbance  code).  An  ad  hoc  three-dimensional 
correction  was  added  to  the  ATRAN2  code  to  intro¬ 
duce  the  cross-flow  effect  along  the  blade.  Cal¬ 
culations  done  with  this  modified  ATRAN2  code  for 
the  rotor-alone  case  showed  the  collapse  of  the 
strong  shock  wave  between  the  azimuthal  positions 
of  150°  and  180°,  and  the  calculations  showei  good 
qualitative  and  fairly  good  quantitative  agreement 
with  experiments.  The  vortex  interaction  flow 
field  for  this  supercritical  condition  was  calcu¬ 
lated  using  this  modified  ATRAN2  code.  The 
results  show  good  agreement  with  experiments  when 
the  vortex  is  near  the  leading  edge  of  the  blade, 
but  seem  to  degrade  progressively  as  the  vortex 
passes  the  airfoil.  In  view  of  the  strong  viscous 
interaction  inherent  in  this  problem,  it  is  not 
surprising  to  see  less-than-perfect  agreement 
between  these  results  and  those  of  experiments. 
Probably  this  is  a  limitation  of  this  approximate 
method.  Nevertheless,  it  ha3  demonstrated  the 
importance  of  the  three-dimensional  effects  at 
these  supercritical  conditions. 


The  computational  procedure  is  very  general 
and  accepts  any  arbitrary  size  and  shape  of  the 
interacting  vortex,  although  the  structure  of  the 
vortex  is  assumed  to  remain  unaltered  by  the 
interaction.  For  comparing  the  present  numerical 
results  with  the  experiments,  it  was  found  that 
the  details  of  the  vortex  structure  were  impor¬ 
tant.  Accordingly,  experimental  data  on  the 
vortex  cylindrical  velocity  distribution  were  used 
to  describe  the  vortex  analytically,  and  this  was 
used  in  the  computations. 

Two  test  cases,  subsonic  and  transonic  condi¬ 
tions  corresponding  to  the  experimental  data,  were 
chosen  for  the  calculations.  The  respective 
blade-tip  Mach  numbers  were  0.6  and  0.8,  and  the 
blade  advance  ratio  was  0.2.  A  comparison  of  the 
numerical  results  for  these  two  conditions  showed 
a  very  distinct  difference  in  the  flow  fields, 
even  in  the  absence  of  the  interacting  vortex. 


A  similar  calculation  done  with  this  approxi¬ 
mate  method  for  the  subcritical  case  showed  almost 
no  three-dimensional  influence,  suggesting  that 
three-dimensional  effects,  if  any,  seem  to  be 
accentuated  in  the  presence  of  shock  waves. 


Appendix 


The  three-dimensional,  transonic  small- 
disturbance  equation  for  a  thin,  high-aspect-rat lo 
rotating  blade,  as  derived  by  Isom, 1:1  is 


A*tt  *  B*xt  -  h  [C1*x  +  C2*x] 


-  9  -  D*  -  E* 

yy  xz  zz 


(At ) 


where 


A  :  M2e2/62/3 
r 


The  results  show  that  for  the  subcritical 
case,  the  unsteady  time-lag  effects  are  negligible 
for  the  rotating  blade  in  the  absence  of  the 
vortex  and  that  the  unsteady  flow  field  with  the 
vortex  Interaction  is  in  very  good  qualitative  and 
quantitative  agreement  with  experiments.  The 
supercritical  case,  however,  is  totally  dominated 
by  strong  (transonic)  shock  waves,  the  consequence 
of  which  is  the  presence  of  strong,  unsteady  time- 
lag  effects,  even  in  the  absence  of  the  vortex 
interaction.  In  addition,  there  are  strong  indi¬ 
cations  of  the  influence  of  three-dimensional 
effects.  The  experimental  data  for  the  rotor- 
alone  case  show  the  collapse  of  the  shock  wave 
between  blade  azimuthal  positions  of  150°  and 
180°.  The  numerical  results,  on  the  other  hand, 
show  the  persistence  of  a  strong  shock  wave  even 
at  180°. 


B  =  2M2c(z  ♦  u'  sin  t)/«2/3 
r 

C,  =  [1  -  M^(z  ♦  V*  sin  t)]/«2/3 

C,  =  1/2( y  ♦  1 )M2(z  ♦  u'  sin  t) 
d  r 

D  r  -2M^eu'  cos  t(z  +  u'  sin  t)/62/3 
E  =  e2/*2/3 

♦  =  #/(Qfi  C«2/3) 

r 

c  =  C/R 

r 

u  *  —  v  /OR 

os  r 

M„  =  nR/a 

r  r  » 


x  =  x/c 
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y  =  y/c 

2  =  5/ft 

r 

t  =  OR  t/R  :  at 
r  r 


In  this  equation,  ♦  is  the  disturbance  potential, 
z  is  the  spanwise  coordinate,  and  x,y  remain 
the  blade  airfoil  coordinates  as  defined  in  Sec¬ 
tion  2  of  the  text.  The  variables  are  scaled  by 
the  airfoil  thickness  ratio  u,  and  by  the  inverse 
of  the  blade  aspect  ratio  e,  following  the  stan¬ 
dard,  small-disturbance  scaling  arguments  for 
transonic  flows21;  namely, 

1  -  M2(z  +  u')  =  O(e)  :0(12/3)  (A2) 

For  the  purposes  of  this  investigation,  the 
reference  conditions  in  Eq.  (A1)  are  defined  for 
the  experimental  rotor  of  Caradonna  et  al.1  as 
follows:  c  is  the  chord  of  the  rotor  blade;  ?B 
is  the  radius  of  the  spanwise  station  where  the 
measurements  were  reported,  that  is,  at 
rB  =  89. 3J  of  the  blade  radius;  Mr  =  (lrg/a2  is 
the  rotational  Mach  number  at  the  measurement 
station;  and  is  the  free-stream  velocity  of 

the  wind  tunnel. 

Spanwise  gradients  that  could  introduce 
strong  three-dimensional  effects  appear  only  in 
the  two  terms,  D*xz  and  £*zz,  in  Eq.  (Al).  An 
approximate,  qualitative  indication  of  how  these 
terms  might  alter  the  two-dimensional  "strip- 
theory"  characteristics  of  the  flow  past  a  rotor 
blade  can  be  obtained  by  considering  a  separation 
of  variables: 

*(x,y,z,t)  =  «(x,y,t)F(z)  (A3) 


Then  Eq.  (Al)  becomes 


Aitt  *  Bixt  -  (C1*x  *  C2*x>x 


V*  =  v*  ♦  qv  (A5) 

where  qy  is  the  prescribed  velocity  field 
induced  by  the  vortex  in  the  absence  of  the  rotor 
blade,  as  explained  in  the  text.  Then  the  govern¬ 
ing  equation  for  ♦  is 

A +  B4  -  (C,«  +  C,*  ($  +  u  )]  -  $ 

tt  xt  lx  2  x  x  v  x  yy 

=  D*V(z)/F(z)  ( A6 ) 

The  left-hand  side  of  Eq.  (A6)  is  unchanged  from 
Refs.  7  and  30,  except  that  the  coefficients  B, 
C^,  and  Cg  now  contain  the  factor  (1  +  u'  sin  t), 
which  was  not  included  before  (see  also 
Ref.  12).  The  value  of  »x  is  known  from  the 
solution  of  Eq.  (A6)  at  the  previous  time-step; 
the  function  F(z)  has  to  be  specified. 

The  boundary  condition  on  the  blade  profile 


ay.  ay. 

*  (y  =  0)  S  ( 1  u'  sin  t)  —  +  e  -rr-  *  j  -  v 


and  the  wake  condition  of  zero  pressure  jump 
across  the  airfoil  branch  cut  is 


[(1  ♦  „•  sin  t)rx  ♦  Ert]wake  =  o 
The  pressure  coefficient  is  given  by 


Cp  =  "  1  -  w’  3in  "  (W  m’  sin  t)2 


Equations  (A6)-(A9)  have  been  solved  using 
the  basic  code  ATRAN2  (Ref.  3d)  with  minor  modifi¬ 
cations,  for  specified  values  of  the  quantity 
F'(z  =  1 )/F(z  =  1). 
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-  i  -  Di  F'/F  -  E*F"(z)/F(z)  (A«) 

yy  x 

where  the  prime  denotes  differentiation  with 
respect  to  z. 

As  a  further  approximation,  we  exclude  the 
tip  region,  thereby  neglecting  the  term  EiF"/F, 
and  assume  that  the  three-dimensional  effects  are 
only  small  perturbations  on  the  basic  two- 
dimensional  solution  at  the  spanwise  station  in 
question,  z  =  z/rB  s  1.  This  allows  us  to  treat 
D*xF'/F  as  an  inhomogeneous  "source"  term,  which 
we  3hall  evaluate  numerically  at  the  previous 
time-3tep,  once  a  suitable  approximation  to 
F'(z)/F(z)  has  been  determined. 

Within  this  framework,  our  previous 
prescribed-disturbance  formulation®’^  for 
vortical  flows  and  the  basic  numerical  algorithm 
for  solving  the  governing  equation  can  be  imple¬ 
mented  with  only  minor  changes.  That  is,  we  write 
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PATH  OF  THE 
MOVING  VORTEX 

Fig.  1  Schematic  of  parallel  blade-vortex  inter-  Fig.  2  Schematic  of  experimental  model  rotor 

action,  in  the  limit  of  A  =  0,  and  definition  of  tests, 

coordinate  system. 
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Fig.  3  Effectiveness  of  the  prescribed-vortex  method  compared  with  the  nonperturbation 
method.  Euler  results  for  airfoil-vortex  interaction:  NACA  64A006  airfoil,  M„  =  0.85, 
•  *  °*f  r  =  0.2,  a0  =  0.05,  and  y0  *  yv  =  -0.26. 
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Fig.  4  Instantaneous  surface-pressure  distributions,  adaptive  grid,  Mach  contours,  and 
pressure  contours  during  airfoil-vortex  interaction  with  a  eonvecting  vortex:  NACA  0012 

airfoil,  M  s  0.8,  a  =  0*,  Re  *  2  *  106  per  meter,  r  *  0.2,  y0  =  yv  *  -0.26. 
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RESULTS 


Fig.  8  Instantaneous  surface-pressure  distributions  for  the  case  of  blade-vortex 


Interaction:  x  0.6,  u  =  0.2,  rB  =  0.893,  r  *  0.243,  yv  =  -0.4. 
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Fig.  9  Instantaneous  surface-pressure  distributions  for  the  rotor-alone  case:  Xj.  =  0.8, 

»  :  0.2,  rB  =  0.893. 
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Fig.  10  Sensitivity  of  the  numerical  calculations  to  the  three-dimensional  correction 
factor  81:  Up  =  0.8,  ti  *  0.2,  and  rB  *  0.893. 
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Fig.  11  Deatons  t  rat  ion  of  the  importance  of  the  three-dimensional  effect  in  the  presence  of 
shock  waves:  Hj.  a  0.8,  it  =  0.2,  and  rB  a  0.893. 
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Fig.  12  Instantaneous  surface-pressure  distributions  without  considering  unsteady  time-lag 

effects:  Mr  a  0.719,  r  a  o.3t,  yv  *  -0.9. 
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Fig.  13  Instantaneous  surface-pressure  distributions  for  the  case  of  blade-vortex  inter 
action,  modified  ATRAN2  solutions:  Hj,  =  0.8,  v  *  0.2,  r  *  0.243,  and  yv  *  -0.4. 


PRESENT  EXPERIMENTS  SURFACE 
RESULTS  (REF.  1) 

-  •  VORTEX 

-  O  UPPER 


Fig.  HI  Typical  results  for  the  case  of  Fig.  13  with  calculations  shifted  by  x„  =  -0.25 
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